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ABSTRACT 

Context. The currently operating space missions, as well as those that will be launched in the near future, (will) deliver high-quality 
data for millions of stellar objects. Since the majority of stellar astrophysical applications still (at least partly) rely on spectroscopic 
data, an efficient tool for the analysis of medium- to high-resolution spectroscopy is needed. 

Aims. We aim at developing an efficient software package for the analysis of medium- to high-resolution spectroscopy of single stars 
and those in binary systems. The major requirements are that the code has a high performance, represents the state-of-the-art analysis 
tool, and provides accurate determinations of atmospheric parameters and chemical compositions for different types of stars. 
Methods. We use the method of atmosphere models and spectrum synthesis, which is one of the most commonly used approaches 
for the analysis of stellar spectra. Our Grid Search in Stellar Parameters (gssp) code makes use of the Message Passing Interface 
(OpenMPI) implementation, which makes it possible to run in parallel mode. The method is first tested on the simulated data and is 
then applied to the spectra of real stellar objects. 

Results. The majority of test runs on the simulated data were successful in the sense that we could recover the initially assumed 
sets of atmospheric parameters. We experimentally find the limits in signal-to-noise ratios of the input spectra, below which the 
final set of parameters gets significantly affected by the noise. Application of the gssp package to the spectra of three Kepler stars, 
KIC 11285625, KIC 6352430, and KIC 4931738, was also largely successful. We found an overall agreement of the final sets of 
the fundamental parameters with the original studies. For KIC 6352430, we found that dependence of the light dilution factor on 
wavelength cannot be ignored, as it has significant impact on the determination of the atmospheric parameters of this binary system. 
Conclusions. The gssp software package is a compilation of three individual program modules suitable for spectrum analysis of single 
stars and individual binary components. The code is highly performant and can be used for spectrum analysis of large samples of stars. 

Key words. Methods: data analysis - Stars: variables: general - Stars: fundamental parameters - Stars: binaries: spectroscopic - 
Stars: individual: (KIC 11285625, KIC 6352430, and KIC 4931738) 


1. Introduction 

Nowadays, spectrum analysis is the main source of precise at¬ 
mospheric parameters and chemical compositions of stars. The 
problem of chemical composition determination in different 
classes of stars is directly linked to the problems of production 
and evolution of chemical elements, stellar and galactic evolu¬ 
tion, etc. Interpretation of the observed atmospheric chemical 
composition in a star gives a general impression about its evo¬ 
lutionary status. For example, it is well known that a star ini¬ 
tially has the chemical composition of the molecular cloud it 
has been formed from. Shortly after the start of nuclear fusion 
in the star, the chemical composition in its central parts under¬ 
goes certain changes. At this stage, the atmospheric composition 
remains quite stable. On the main-sequence, certain processes 
may occur in the star that will cause the exchange of matter be¬ 
tween the stellar interior and the atmosphere. Examples are the 
products o f the CNO-cyc le in the atmosph eres of massive O-B 
stars (e. g.. iMaederll 1 983alfhL lMassevll2003l) . or lithium depletion 
in the atmospheres of F-K dwarf stars (e.g., ISoderbloml2015t 
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lAndrassv & Snruitll2()15f iDelgado Mena et all 120151) . Some of 
the processes that can cause the above mentioned mixing of the 
material are convection or turbulent diffusion. 


At early stages of giant evolution, the stars generally go 
through the phase of deep convective mixing which has cer¬ 
tain impact on the atmospheric composition of light (and later 
on heavier) chemical elements. Detailed studies of atmospheric 
chemical composition can also reveal signatu res of mass loss 
in massive stars (e.g., Chiosi & Maeder 1986) and episodes o f 
mass transfer in close binary systems (e.g.. IBerdvuginalll994l) . 
Chemically peculiar stars are another example of how detailed 
studies of atmospheric chemical composition allow for a better 
understanding of different physical processes occurring in stars 
(e.g JPohnl et al]l200l iKochukhov & BagiiuiollIOOhh . 

Tha nks to the launche s of seve ral space missions su ch as 
MOST (Walk er et al.ll2003l). C oRoT (lAuvergne et al.ll2009l) . and 
Kepler (Gilliland et al. 2010), the astrophysical fields such as, 
e.g., asteroseismology and exoplanetary science has been rev¬ 
olutionized during the past decade. High quality photometric 
data obtained from space revealed countless number of interest¬ 
ing physical effects in different types of stellar objects and led 
to very interesting discoveries. Needless to say, however, that 
both fields still heavily depend on high-quality, high-resolution 
spectroscopic data, and a significant fraction of the analyses 
rely on the interpretation of combined space-based photomet- 
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Table 1 . Stellar atmosphere models computed with the LLmodels code 
for £ = 2 km s -1 . The table is re-produced from lTkachenko et alj d2012l) 


Parameter, step width 

[M/H] 

A [M/H] 

Teff(K) 

ATeff(K) 

logg 

Alogg 

-0.8-+0.8 

0.1 

4500-10000 

100 

2.5-5.0 

0.1 

10 000-22 000 

250 

3.0-5.0 

Total number of models: 



41888 


ric and ground-based spectroscopic data. The recently launched 
Gaia mission dPerrvman et al.ll200ll) and missions such as TESS 
(iRicker et al.ll2014b and PLATQ2.0 dRauer et al.ll2014l) that will 
be launched in the near future, will provide high-quality data for 
a few million of stellar objects suitable for spectroscopic follow¬ 
up from the ground. As such, efficient and fast tools are particu¬ 
larly needed for the analysis of ground-based spectroscopic data 
of single and multiple stellar objects. 

Nowadays, the method of atmosphere models and spec¬ 
trum synthesis is a dominant approach in the analysis of high- 
resolution spectra of single stars and those in binary systems. 
One of the advantages of the spectrum synthesis over the tradi¬ 
tional methods that rely on the calculation of equivalent widths 
is that the effects of line blending can be accurately taken into 
account. This in turn minimizes the uncertainties in the funda¬ 
mental stellar parameters and chemical abundances determina¬ 
tion. One example of widely used software for spectrum anal¬ 
ysis of single stars based on the spe ctrum synthesis method is 
the Spectroscopy Made Easy (SME, IValenti & Piskuno~vl[1996i) 
package. 

In this paper, we present a new Grid Search in Stellar 
Parameters (gssp) software package for spectrum analysis of 
high-resolution spectra of single stars and those in binary sys¬ 
tem. In the next sections, we describe the implemented method¬ 
ology and test the code both on simulated and real stellar spectra. 

2. Methodology 

In this section we discuss in all necessary detail the methodol¬ 
ogy implemented in the gssp software package. This includes the 
description of each of the three program modules and the results 
of the test runs. In summary, the gssp_single module is designed 
for the analysis of single star spectra and the disentangled light 
diluted spectra of the individual components of double-lined 
spectroscopic binary system, assuming wavelength-independent 
light dilution factor. In this case, the individual binary compo¬ 
nents are interpreted independent of each other, just as they were 
single stars, unlike the method implemented in the gssp_binary 
module which is specifically designed for simultaneous inter¬ 
pretation of the binary components’ disentangled spectra, by 
taking into account wavelength dependence of their light di¬ 
lution factors. A similar methodology is implemented in the 
gssp_composite module, where the fitting of the disentangled 
spectra is replaced by the fitting of the observed composite spec¬ 
tra, and the wavelength dependence of light dilution is also as¬ 
sumed. 

2.1. Basic methodology 

Although there are small differences in the realization of the in¬ 
dividual algorithms, the general methodology is the same for all 
three modules. The software package is based on a grid search in 
the fundamental atmospheric parameters and (optionally) indi¬ 


vidual chemical abundances of the star (or binary stellar compo¬ 
nents) in question. We use the method of atmosphere models and 
spectrum synthesis, which assumes a comparison of the observa¬ 
tions with each theoretical spectrum from the grid. For calcula¬ 
tion of synthet ic spectra, we use the S ynthV LTE-based radiative 
transfer code (lTsvmballll996l) and a grid of atmosphere mod - 
els precomputed with the LLmodels code (IShulvak et al .1 12004). 
Our grid of models covers large ranges in all fundamental at¬ 
mospheric parameters and is freely distributed together with the 
software package itself. The summary of the available models 
is given in Table Q] In fact, the gssp package is compatible with 
any kind of atmosphere model grid as long as the models are 
provided in the Kurucz format. The authors also possess a grid 
of Kurucz rnodelf] (tKuruczll 19931) which have been interpolated 
in all fundamental parameters to match the resolution of our 
LLmodels grid. The interpolated grid of Kurucz models is avail¬ 
able upon request. 

We allow for optimization of five stellar parameters at a 
time: effective temperature T e ff, surface gravity log g , metallic- 
ity [M/H], microturbulent velocity £■', and projected rotational 
velocity v sin i of the star. The synthetic spectra can be com¬ 
puted in any number of wavelength ranges, and each consid¬ 
ered spectral interval can be from a few angstroem up to a few 
thousand angstroem wide. As long as the global metallicity of 
the star is determined/known, the [M/H] parameter can be re¬ 
placed in the grid by the abundance of an arbitrary chemical el¬ 
ement. The individual abundances have to be iterated element 
by element, thus there is no option to optimize abundances of 
more than one element at the same time. From our experience, 
any reasonable deviations of the individual abundances from the 
global atmospheric metallicity (within about 0.5 dex) have little 
(within typical lcr error bars) to no influence on abundances of 
other chemical elements. An exception needs to be made how¬ 
ever for the elements showing a large amount of lines in their 
spectrum: should such an element be found to show significant 
over/underabundance, all chemical element abundances iterated 
before must be re-determined, taking into account the fact that 
one of the dominant chemical elements in the spectrum shows 
peculiarities. For this reason, we advice the user to start with 
chemical elements having the largest amount of lines in the spec¬ 
trum when determining the detailed chemical composition of the 
star in question. 

The grid of theoretical spectra is built from all possible com¬ 
binations of the above mentioned parameters. Each spectrum 
from the grid is compared to a priori normalized observed spec¬ 
trum of the star and the x~ merit function is used to judge the 
goodness of fit. The code delivers the set of best fit parame¬ 
ters, the corresponding synthetic spectrum, and the ASCII file 
containing the individual parameter values for all grid points 
and the corresponding x 2 values. The x 1 we deliver is the re¬ 
duced^ 2 , that is the^ 2 value normalized to the number of pixels 
across the observed spectrum minus the number of free parame¬ 
ters. Based on the^ 2 statistics and assuming normal distribution, 
we also compute 1 cr uncertainty level in terms of x 2 and output 
this value along with the x 2 distributions. We also account for 
possible global-scale imperfection in the normalization of the 
observed spectrum by means of a scaling factor that is computed 
from the least-squares fit of the synthetic spectrum to the obser¬ 
vations and is applied to the latter. The value of the scaling factor 
is provided in the final x 2 table along with the other grid search 
parameters. 


1 http://kurucz.harvard.edu/grids.html 
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Fig. 1. Diagram illustrating the algorithm implemented in the 
gssp.single software module. The two identical branches showed in the 
diagram indicate multiprocessing. See text for more details. 


For the analysis of binary stars, whether the disentangled 
spectra or the observed composite spectra are used for the char¬ 
acterization of the system, the percentage contribution of a stel¬ 
lar component to the total light of the system becomes one of the 
key parameters. This light contribution is often referred to as a 
light dilution factor (designated as f throughout this paper) and 
effects the depths of lines in the spectra of both binary compo¬ 
nents. This correction factor for the line depths has to be taken 
into account in the analysis, and is ideally determined along with 
the other atmospheric parameters of the star. There are two basic 
ways of accounting for the light dilution factor in the spectro¬ 
scopic analysis: i) assuming a constant scaling of the line depths 
across the entire wavelength range (unconstrained fitting of the 
disentangled spectra), and ii) taking into account wavelength- 
dependence of the light dilution due to possible difference in 
spectral type and/or luminosity class for the two stars forming 
a binary_system (constrained fitting of the disentangled spectra, 
iTamaio et al.11201 l~l) . Which of the two methods to use depends 
on the analysis approach, on whether the individual binary com¬ 
ponents are characterized independent of each other or simul¬ 
taneously by means of fitting their disentangled spectra. More 
details on both methods are given in Sect. 2.2 and 2.3. 

The grid search method is generally more CPU time de¬ 
manding than any of the optimization algorithms but has a big 
advantage that it guarantees that the global minimum solution 
will be found as long as the considered parameter range is large 



Fig. 2. Quality of the fit of the simulated data of KIC 11285625 (model 
lb in Table lA.it with the gssp.single software package. The light dilu¬ 
tion factor is 0.6/0.4 for the primary/secondary component. The simu¬ 
lated observed spectrum is shown with light grey lines, the black solid 
line refers to the best fit synthetic spectrum. The spectrum of the sec¬ 
ondary was vertically shifted by a constant value for better visualiza¬ 
tion. 


enough. We use the OpenMPQ distribution to parallelize our 
code, which solves the CPU time-related problem and makes the 
code fast and highly performant. At least for one of the software 
modules, there seems to be a linear relation between the num¬ 
ber of used CPUs and the calculation time. However, there still 
exists an upper limit for the number of CPUs to be used. The 
limit strongly depends on the configuration of the system and 
is likely to occur when the system load becomes large enough 
to start effecting the performance of the code. We will provide 
rough estimates of the calculation times below, when discussing 
each of the gssp package modules individually. 


2.2. gssp_single module 

In this particular module we treat the observed spectrum as the 
one of a single star. The module is thus applicable to spectra of 
single stellar objects as well as to the disentangled spectra of 
multiple stellar systems. In the latter case, we assume that the 
light of the star in question is diluted by a certain wavelength- 
independent factor, and the disentangled spectrum of a star is 
represented as follows 


Rf s = 


N -1 

f v j + £ f; 

j= i ’ 


£ F 

k=l 


c 

v,k 



(i) 


where I v is the specific intensity of the star at a given wave¬ 
length/frequency, X is the radius of the star, and d is the distance 
to the object. After a few trivial mathematical operations, for ei¬ 
ther of the components of the binary system we obtain in terms 
of line depths 


„dis 


r* 


1 + a 


iff' 


1 v,comp 


x com 


2 http://www.open-mpi.org/ 
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Table 2. Results of the spectrum analysis of both stellar com¬ 
ponents of KIC 11285625 with the gssp_single and gssp_binary 
softw are packages (s e e text for details). The parameters derived 
bv iDebosscher et all f 2013 1) are given as a reference; the log g 
values and the radii radio were derived by the authors from light 
curve solution. 


Parameter 

Unit 

GSSP_SINGLE 

primary secondary 

GSSP_B INARY 

primary secondary 

T e s 

K 

7 200/®° 7 300/;°° 

6980/“ 7100/!™ 

log 8 

dex 

3.91/°“} 4.18 (fixed) 

3.82/°;! 4.18 (fixed) 


km s _1 

1 . 19 /°;° 0 . 0 + °- 6 

O 

K> 
l + 
© P 

O 

O 

+ 

© 

0 

v sin i 

km s~* 

14-5/1“ 8 . 6 / 1 ° 

15.6/!;? 8.9/!;? 

[M/H] 

dex 

0 9C+O.O8 A AI+O.IO 

-0.08 _W ^ A -0.10 

_0 43+ 0 - 09 _0 37 +°- 15 
-0.09 ^ ,J/ -0.15 

R 1 /R 2 (f) 


a 54+0.05 a 27+0.02 
u,u -0.05 

1 62 +016 

1-U -0.15 


Parameters from Debosscher et al. (2013) 



primary 

secondary 

7eff 

K 

6 960+ 100 
u ^ UVA -100 

7 195 +200 

logg 

dex 

3.97 (fixed) 

4.18 (fixed) 


kms -1 

0 95 +0 - 30 

-0.30 

0 09 +0 - 25 

-0.25 

v sin i 

km s -1 

14 2 +1 - 5 

A ^‘ —1.5 

8 4 +L5 
—1.5 

[M/H] 

dex 

-0 49 +0 - 15 
-0.15 

-0 37 +0 - 30 
u,J/ -0.30 


R,/R 2 1.441°;“ 

Note: dilution factors f are reported for the gssp_single mode 


Table 3. Same as 'fable [2] but for KIC 6352430. The reference 
parameters are those from lPapics et ali (i2013l) . 


Parameter 

Unit 

GSSP_S INGLE 

primary secondary 

GSSP_BINARY 

primary secondary 

7rfr 

K 

13 200 +10 ° 7 065 +10 ° 

LJ 4UU_ 100 / UUJ_ 100 

13 060/!!? 7 340/?“ 

logg 

dex 

4 O 7 +0 04 3 70+0.20 

^•^-0.04 ' -0.20 

4 03 +0 05 3 60 +0 - 56 

-0.05 -0.56 


kms -1 

1 1 a+0.53 0 '7C+0.50 

1,1 -0.58 7 -0.50 

1 9Q+0.43 4 5Q+0.90 

1 , yu - 0.41 -0.90 

v sin i 

kms -1 

11 Q+4.1 Q c+1.0 

-3.8 ® -0.9 

71.4/?;? 7.0/?;? 

[M/H] 

dex 

-0.08/2! -0-24/gJg 

-0.19/°! 0.01/?;!® 

R 1 /R 2 ( f ) 


0-95/°;°! 0.06/°;“ 

1 . 86 /?;!! 


Parameters from Panics et al. (2013) 



primary 

secondary 

7eff 

K 

12 810/?“ 

6 805/!?? 

logg 

dex 

4.05/™? 

4.26/?;!? 

? 

km s -1 

2.0 (fixed) 

2.0 (fixed) 

v sin i 

km s -1 

69-8/?;° 

9 8 +1 ° 

[M/H] 

dex 

-0.13/2;“ 

-0 33 +01 ° 

-0.10 


Note: dilution factors f are reported for the gssp_single mode 


where superscripts “th” and “c” refer to the theoretical spectrum 
and continuum intensity, respectively, and subscript “comp” 
points to the companion star. The factor fj = 1 /(1 + a) is the 
light dilution factor briefly outlined in Sect. 2.1. In this particular 
case of interpreting the individual binary components’ spectra 
independent of each other, nothing is known about the contin¬ 
uum intensity ratio of the two stars and the light dilution factor 
f is assumed to be wavelength independent. In practice, each 
synthetic spectrum from the computed grid is represented in the 
form of Eq. [2] and compared to the observed disentangled spec¬ 
trum on the scale of the latter. 


Table 4. Same as Table[2]but for KIC4931738. The reference 
parameters are those from lPapics et at . ( 201 3). 


Parameter 

Unit 

GSSP_SINGLE 

primary secondary 

GSSP_BINARY 

primary secondary 

7eff 

K 

14120/??? 11975/??? 

14150/??? 11780/1? 

logg 

dex 

4 04+ 010 4 40+ 011 

^•^-0.09 ^ u - 0.11 

415 + 0.07 450 +O.I 6 

^• A -0.09 -0.16 

£ 

km s“ 

0.40/?;?? 0.88/?;?? 

0.00 +ll ° 0.55 +l - 45 

v sin i 

km s -1 

11.3/!;? 6.5/JI 

H.2/?;? 6.6/[;? 

[M/H] 

dex 

-0 15 +009 0 06 +009 

U-A - 0.10 - 0.10 

-0.19/?;?! 0.14/?;!? 

R 1 /R 2 (I) 


0 79 + 0.03 0 04 + 0.02 

w,/ -0.03 -0.02 

1-26/?;?? 


Parameters from Panics et al. (2013) 



primary 

secondary 

T e s 

K 

13 730+ 200 

LJ /2JU_ 20 q 

11 370 +250 

11 -250 

logg 

dex 

0 Q 7 +O.O 5 
' -0.05 

4 / 27 + 0.10 
^• J/ -o.io 

£ 

km s _1 

2.0 (fixed) 

2.0 (fixed) 

v sin i 

km s _1 

10 5 +1 ° 

AU - 1.0 

6-6!| : ? 

[M/H] 

dex 

-0 24 +010 
u> ^-0.10 

+ 0 . 10 /?;!? 


Note: dilution factors f are reported for the gssp.single mode 


Figure Q] illustrates the multiprocess algorithm implemented 
in the gssp_single software module. In the first step, the code 
reads in all necessary information provided by the user in the 
configuration file. We provide an example of the input file with 
detailed description of each entry in the Appendix. In the next 
step, the code sets up the grid and checks whether the required 
grid of atmosphere models is available. Should any of the models 
not exist, the user will be notified about the need to revise the in¬ 
put set up, and the detailed report is saved in one of the log-files. 
At this point, the observed spectrum can be optionally cross- 
correlated with the first synthetic spectrum from the grid, and 
the radial velocity (RV) is computed as the first order moment of 
the cross-correlation function. This way, the code accounts for a 
possible RV shift of the observed spectrum with respect to the 
laboratory wavelength, should the user request that. As soon as 
the grid has been set up (depending on the grid size, this might 
take up to half a minute of time), the code enters a multipro¬ 
cessing mode and performs the actual calculations. The SynthV 
code is executed on each of the available CPUs and provides a 
synthetic spectrum for a given set of atmospheric parameters and 
in the requested wavelength range. The code computes specific 
line and continuum intensities for different positions of the stel¬ 
lar disk, thus providing an accurate treatment of the limb dark¬ 
ening effect. As soon as the spectrum is computed, a convolution 
program is executed to perform the disk integration and a con¬ 
volution of the spectrum with projected rotational and macrotur- 
bulent velocities, and the resolving power of the instrument. The 
convolution code outputs the normalized synthetic spectrum as 
well as the line and continuum fluxes. The above described cal¬ 
culations may take from a few seconds up to a few minutes, de¬ 
pending on the spectral type and considered wavelength range. 
Once the convolved synthetic spectrum is released, there are two 
options for the code to contiue with (see Fig.|T}: immediate com¬ 
parison with the observations (the case of the spectrum of a sin¬ 
gle star), or dilution according to Eq.[2]and then the comparison 
with the observations (the case of the disentangled spectrum of 
a binary component). In either case, the final output of an indi¬ 
vidual process is the goodness of fit of the observed spectrum 
in terms of the x 2 merit function value, whereas the computed 
synthetic spectrum is immediately deleted. As soon as the cal- 
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Fig. 3. Comparison between the observed disentangled spectra (light 
grey lines) of both stellar components of the KIC 4931738 system and 
the best fit synthetic spectra (black lines) computed with the gssp_single 
software module. The spectrum of the secondary was vertically shifted 
by a constant value for better visualization. 


culations are finished on all the CPUs, the main process merges 
individual^- 2 tables, finds the minimum in y 1 , computes the best 
fit synthetic spectrum, and outputs the set of best fit parameters. 
We also provide the^ 2 value corresponding to the 1 <r uncertainty 
level computed from x 2 statistics and assuming normal distribu¬ 
tion. This way, the projection of all data points on the parame¬ 
ter in question can be computed afterwards from the x 2 “master 
file”, and the lcr uncertainties are obtained in a straightforward 
manner by fitting the x 2 distribution with a polynomial function 
and searching for its intersection points with the lcr level in x 1 
provided by the gssp_single module. 

The gssp_single module is reasonably fast and the final so¬ 
lution is obtained within 5-6 minutes by running the code on 8 
CPUs, and assuming a typical grid size of 2 000-3 000 synthetic 
spectra and a 1 000 A wide wavelength range. The calculations 
may take a bit longer when fitting disentangled spectra, as there 
is an additional parameter (light dilution factor) to account for. 
The extra calculation time typically does not exceed 10% of the 
computation time in a single star spectrum mode, unless the grid 
resolution in the dilution factor is made unreasonably high. 

The gssp_single software module is in fact_ an exten¬ 
sion of_the_original_version of the gssp code dLehmann et alJ 
1201 It iTkachenko et al.l l2012h to make it compatible with the 
disentangled diluted spectra of stellar components in binary 
systems. The module was tested on simulated data of sin¬ 
gle stars]_Jn_addition, we refe^_the_^eader_tothe_papers 
of iLehmann et all ( 201 ll), Tkache nko et al. (1201 2l 201 3alfbt). 


iLam pens et al.1 d2013l) . IPanics et al] (12014 ). I Van Reeth et al J 
d2015l) for several examples of the code application to the spectra 
of single stars, including some well studied objects like Vega. 

For tests on the light diluted spectra of individual binary 
components, we decided not to simulate the spectra of arbitrary 
objects but instead used the binary com ponents of KIC 6352430 
dPdnics et al .1120131) and KIC 11285625 dDebosscher et al.ll2013l) 


as prototype stars. In our simulations, we assumed different dilu¬ 
tion factors as well as different values of S/N ratio; the summary 
of the obtained results is given in Table I A. II Overall, we could 
recover the parameters of the input spectra assumed in all our 
simulations; the largest deviations, as well as lcr uncertainties, 


have been observed for low S/N ratio spectra. The quality of the 
fit to one pair of spectra is illustrated in Fig. [2] The only set up 
where we have encountered certain difficulties in recovering the 
assumed fundamental parameters refers to the case of peculiar 
stars (see Table lAAl the column indicated as “gssp_single”). In 
the simulations of KIC 11285625, we assumed an over/under¬ 
abundance of iron/silicon by 0.3/0.5 dex with respect to the 
global metallicity of the star for the primary/secondary com¬ 
ponent. For the KIC 6352430 system, the primary star was as¬ 
sumed to be peculiar in helium (+0.2 dex compared to the solar 
composition), whereas the secondary component was set up as 
a normal star with its own global metallicity. Since in the un¬ 
constrained fitting the spectra of individual components are con¬ 
sidered separately, there is no way the results obtained for one 
of the components will influence the parameters of the compan¬ 
ion. For this reason, we give no parameters for the secondary of 
KIC 6352430 in Table IA.4I as they will be the same as the ones 
depicted in the second column of Table 0 We find that when a 
chemical element dominating the observed spectrum of the star 
shows significant over-/underabundance (iron in the primary of 
KIC 11285625), the global metallicity value gets affected and 
impacts several other fundamental parameters as well (e.g., mi- 
croturbulent velocity and surface gravity). Although this is not 
an unexpected result and the majority of the parameter values 
are still within lcr uncertainties from the assumed values, the 
fact calls for a certain attention. 

After testing the code on simulated data, we have applied 
the gssp_single module to the real disentangled spectra of both 
stellar components of the KIC 11285625, KIC 6352430, and 
KIC 4931738 systems. The latter object was also analysed by 
iPanics et al.l ( 2013! ) and consists of two stars of similar spectral 
types, luminosity classes, and rotation rates. The results of our 
analyses for all three binary systems are depicted in Tables [2] 
0 and [4] respectively, in the columns indicated as “gssp_single”. 
The parameters obtained by us are slightly different from those 
presented in the original studies, though the majority of the pa¬ 
rameters still agree within the reported lcr uncertainties. The 
obtained discrepancies for KIC 6352430 and KIC 4931738 may 
be associated with the fact that the microturbulent velocity was 
fixed to 2 kms -1 in the original study bv IPanics et al .l (UOI 3t ). 
whereas it was set as a free parameter in our analysis. In ad¬ 
dition, Panics et~akl (1201 3l) assumed some dilution factor values 
coming from the preliminary fit of several observed composite 
spectral lines and used those to re-normalize the disentangled 
spectra of both binary components to the individual continua. In 
our case, both factors were set as free parameters and optimized 
along with other fundamental parameters for each of the compo¬ 
nents. As an example, in Figure 0 we show the quality of the fit 
of the observed disentangled spectra of both components of the 
KIC 4931738 system by the best fit synthetic spectra computed 
from the parameters depicted in the first two columns of Table0 

As it has been mentioned bv lDebosscher et all (120131) in the 
original study, the spectroscopic data of KIC 11285625 have low 
S/N, which naturally propagates to the quality of the disentan¬ 
gled spectra. The secondary component, which contribution to 
the total light amounts to ~30%, suffered the most and its disen¬ 
tangled spectrum has rather low S/N and shows pronounced un¬ 
dulations in the local continuum. We anyway attempted to fit the 
spectrum of the secondary setting T e g, log g, / ; , v sin i, and [M/H] 
as free parameters in our analysis. We ended up with a clearly 
unreliable solution, in particular the value of log g was found to 
exceed 5.0 dex, which is not expected for the main-sequence star 
of spectral type F. Thus, following the p rocedure adopted in the 
original study of lDebosscher et all ( 2013 ). we fixed surface grav- 
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Fig. 4. Diagram illustrating the algorithm implemented in the 
gssp_binary software module. The two identical branches showed in the 
diagram indicate multiprocessing. See text for more details. 


ity of the secondary to the value obtained from the light curve 
solution, and optimized the four remaining fundamental param¬ 
eters. Since the disentangled spectrum of the primary has higher 
S/N and suffered less from artifacts of spectral disentangling, we 
optimized all five fundamental parameters for this star, with the 
particular goal to compare the spectroscopic log g value with the 
one obtained from the light curve solution. The results of the 
analysis are summarized in Table [2] and show good agreement 
between the spectroscopic and photometric values of log g for 
the primary component, and overall agreement for all parame¬ 
ters of both binary components with the original study. 

2.3. gssp_binary module 

This module has been specifically designed for fitting disentan¬ 
gled spectra of both binary components simultaneously. The pro¬ 
cedure is in a sense similar to the one of “constrained fitting” 
suggested by iTamaio et ab (1201 li) . were one assumes that the 
sum of the two light factors is identical to unity, and thus the 
change in light dilution of one of the components affects the 
amount of diluted light for the companion star. Instead of assum¬ 
ing wavelength-independent light dilution for each of the com¬ 
ponents, we optimize the ratio of the radii of the binary compo¬ 
nents which is obviously the same for both stars. Starting from 
the definition of disentangled spectra of individual binary com¬ 
ponents given in Eq. |T| and dividing all terms on the right hand 


side by the product between the continuum intensity of the sec¬ 
ondary and its squared radius, we obtain for the primary and 
secondary components of a binary 
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Just as in the formulation of iTamaio et al.l (1201 It) , the sum of the 
two dilution factors is identical to unity. By optimizing the ratio 
of the radii, we take into account the wavelength-dependence of 
the light dilution factor through the ratio of continuum intensi¬ 
ties of the two stars. Obviously, any change in the atmospheric 
parameters of one of the binary components will influence its 
continuum intensity, which will in turn impact the light dilution 
factor for the companion star. In practice, disentangled spectra of 
two binary components are analyzed simultaneously by scaling 
synthetic spectra from the corresponding grids and comparing 
them to the observations on the scale of the latter. 

Figure [4] visualizes the algorithm implemented in the 
gssp_binary software module. Although the general principle is 
the same as for the gssp_single module, there are a number of 
differences that we discuss here in more detail. After reading all 
necessary information from the input file and checking the at¬ 
mosphere models for availability, the code sets up the grid of 
synthetic spectra that need to be computed. As we have to deal 
with two components simultaneously, the algorithm of setting up 
the grid is different from the one implemented in the gssp_single 
module. To supply the code with maximum efficiency, we first 
build individual grids for both components of a binary system, 
and in the second step, a unique grid common for both stars is 
build by excluding any possible overlap in the parameters be¬ 
tween the individual grids. This step is needed to avoid unneces¬ 
sary repetitive calculations of theoretical spectra. For example, 
in the case of a binary system with two identical stars and thus 
identical initial parameter grids, only the primary’s grid will be 
computed and used for both binary components afterwards. The 
individual grids are used at this step to build all possible com¬ 
binations of primary and secondary spectra, which will be pro¬ 
cessed at a later stage when looking for the best fit solution. For 
example, three grid points for each of the binary components will 
result in nine combinations to be processed in total, regardless 
of any possible overlap between the individual grids. As soon 
as the grids were set up (takes less than a minute of time), the 
code enters the parallel mode and starts the actual calculations 
of synthetic spectra as described in Sect. 2.1. 

The normalized synthetic spectra as well as the correspond¬ 
ing continuum fluxes are stored in a database, which is an un¬ 
formatted direct access Fortran 90 file. The advantage of using 
such a file over storing the data in an ASCII file is that any I/O 
operation with such a file is much faster as there is no conver¬ 
sion from machine code to a readable format. By storing the 
data in such a file we also avoid any possible problems with in¬ 
sufficient RAM memory, which might occur when keeping the 
(sometimes very large) grid of synthetic spectra in the memory. 
As soon as the database of spectra is produced, it is first made 
available to all CPUs involved in the calculations, and then the 
code enters the parallel mode for the second time. At this stage, 
every pre-build combination of primary and secondary spectra is 
processed, which involves dilution of both synthetic spectra ex¬ 
tracted from the database and comparison with the observations 
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on the scale of the latter. After processing all possible combina¬ 
tions, they' 2 files from individual CPUs are merged into a “mas¬ 
ter file”, and the best fit parameters and synthetic spectra for 
both binary components are released. Similar to the gssp_single 
software module, the lcr level in x 1 is provided so that the cor¬ 
responding uncertainties in all the parameters can be computed 
afterwards. 

Despite a very large number of combinations that have to be 
processed, the gssp_binary module is reasonably fast. Assuming 
the same set up as the one described in Sect 2.1-8 CPUs and 
a 1 000 A wide wavelength range, - and ~500 synthetic spectra 
in a grid, the final solution is obtained within 15-20 minutes. To 
provide an estimate on the number of combinations, the above 
mentioned grid of ~500 spectra coupled with only three grid 
points in the radii ratio results in about 200000 combinations 
to be processed by the code. 

The gssp_binary software module has been extensively tested 
by us both on simulated and real spectra of double-lined spec¬ 
troscopic binary stars. Similar to the tests described in Sect 2.1, 
KIC 11285625 and KIC 6352430 were used as prototype stars 
for our simulations, and the observed disentangled spectra of the 
KIC 6352430, KIC 11285625, and KIC 4931738 systems were 
used to test the code on data of real binary systems. The results 
of the code application to the simulated data are summarized in 
Table IA.2I Similar to the tests performed in Sect. 2.1, we also 
applied the gssp_binary module to the simulated spectra show¬ 
ing anomalies in the abundances of certain chemical elements. 
The results of this application are summarized in Table IA.4I in 
the columns indicated as “gssp_binary”. The conclusions from 
our test runs on the simulated data are essentially the same as in 
Sect. 2.1: i) the input parameters are recovered in the majority of 
the cases, ii) the largest deviations from the input values as well 
as the largest uncertainties are observed for low S/N data, and iii) 
anomalies in the abundances of chemical elements dominating 
the observed spectrum significantly affect the determination of 
the global metallicity, which in turn triggers deviations in other 
fundamental parameters of the star. Since in the constrained fit¬ 
ting mode spectra of both binary components are fitted simulta¬ 
neously, any changes in the parameters of one of the components 
will likely affect the parameter determination for the companion 
star, due to changes in the light dilution factor. The quality of the 
fit of the simulated data is essentially the same as the one shown 
in Fig. [2] for which reason it is not illustrated here. 

The results of the application of the gssp_binary software 
module to the real spectra of KIC 11285625, KIC 6352430, and 
KIC 4931738 are summarized in Tables [2] [3] and |4j respec¬ 
tively (column indicated as “gssp_binary”). Significant devia¬ 
tions are observed for some of the parameters from those ob¬ 
tained with the gssp_single module (unconstrained fitting) for 
the KIC 11285625 system. This is likely explained by the insuf¬ 
ficient quality of the data for the secondary component, whose 
parameter determination in turn impacts the solution for the pri¬ 
mary component. Indeed, the results of the unconstrained fit¬ 
ting with the gssp_single module (first two columns in Table [2} 
suggest the sum of the dilution factors to be well below unity. 
According to the definition of disentangled spectra from Eq. [3] 
this sum is identical to unity in the constrained fitting mode, thus 
the gssp_binary code is in a sense forced to look for an alternative 
solution that would satisfy the above mentioned criterion. The 
change in the light dilution of each of the components in turn 
triggers the differences in fundamental parameters compared to 
the gssp_single solution. 

Certain deviations in the metallicity and effective tempera¬ 
ture of the secondary are observed for the KIC 6352430 when 



Fig. 5. Ratio of the continuum intensity of the primary to the one of the 
secondary in the KIC 6352430 system. The models were computed for 
the fundamental parameters listed in the last two columns of Table [3] 


compared to the solution obtained with the gssp_single module. 
Given that the latter suggests unity for the sum of the light fac¬ 
tors, the observed discrepancies are likely due to wavelength de¬ 
pendence of the light dilution factors. Indeed, the primary com¬ 
ponent is twice as hot as its companion star, suggesting that 
the continuum intensity ratio might vary significantly within the 
considered 1 000 A wide wavelength range. Figure 0 illustrates 
the variation of the ratio of the continuum intensity of the pri¬ 
mary to the one of the secondary as a function of wavelength, 
in the range between 4 750 and 5 750 A. The ratio changes from 
~4.9 at AA 4750 A to ~3.9 at AA 5 750 A. For comparison, the 
variation in the same quantity for the KIC 11285625 system con¬ 
sisting of two similar stars (see Table O) does not exceed 0.2% 
of the value of 0.955 measured at 5 250 A. 

2.4. gssp_composite module 

The gssp_composite software module was specifically designed 
for fitting composite spectra of double-lined spectroscopic bi¬ 
nary systems. We refer the reader to Fig. |4] and the description 
there for details on the implemented algorithm as it is very simi¬ 
lar to the one implemented in the gssp_binary module. The only 
two differences in the gssp_composite algorithm are: i) there is a 
possibility to set the radial velocities of individual binary com¬ 
ponents as free parameters, and ii) instead of comparing theo¬ 
retical spectra to the disentangled observed spectra of both com¬ 
ponents, all possible combinations of synthetic primary and sec¬ 
ondary spectra from the computed grid are used to build com¬ 
posite theoretical spectra of a binary which are then compared 
to the observed spectrum on the scale of the latter. Composite 
theoretical spectra are represented in the form 

/v.i^T + AAX? aRf+Rf 

Rcomp — yy 1 . VJ 2 — TT ’ (5) 

I , zvr + / 'a~ 1 +a 

v,l 1 v,2 2 

where a is defined the same way as in Eq. [3] 

We have tested the gssp_composite module on simulated 
composite spectra of binary stars, where the KIC 11285625 and 
KIC 6352430 binaries were again used as the prototype systems. 
The results of our test runs are summarized in Table lA.31 quality 
of the fit to one of the spectra is illustrated in Fig. [6] As for the 
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Wavelength (Angstroem) 


Fig. 6. Quality of the fit of the simulated composite spectrum of 
KIC 11285625 (model la in Table liOl l with the gssp_composite soft¬ 
ware package. The radii ratio is 1.45, in agreement with the value re¬ 
ported bv iDebosscher et ah ( 2013 1. The simulated observed spectrum 
is shown with a light grey line, the black solid line refers to the best fit 
synthetic spectrum, the individual spectra of the primary/secondary are 
shown by the red dashed/blue dotted lines. The individual, RV-shifted 
spectra of the binary components were vertically shifted for better visu¬ 
alization. 


two previously described software packages, we could recover 
the input parameters for the majority of the models. However, 
we find that a combination of low S/N and a small light con¬ 
tribution from either of the binary components becomes a real 
bottleneck for fitting the composite spectrum of a binary sys¬ 
tem. The only reasonable estimates in this case are obtained for 
the radial velocities of stars, provided both components show 
a sufficient amount of metal lines in their spectra and the lines 
are well separated in velocity space. For the rest of the parame¬ 
ters fitting the composite spectra fails, for which reason there are 
no entries for the secondary components of KIC 11285625 and 
KIC 6352430 in Table IA.3I at S/N of 60 and large values of the 
radii ratio. An assumption that one or both stellar components 
of a binary show significant anomalies in their spectra leads to 
the same results than have been obtained with the gssp_binary 
software package (see Sect. 2.2). For a limited amount of mod¬ 
els, we also performed test runs assuming that spectral lines of 
individual spectral components overlap in the velocity space. 
To our surprise, even in this case all initially assumed funda¬ 
mental parameters could be reasonably recovered, provided suf¬ 
ficiently high S/N (>150) of the input data. No application of 
the gssp_composite module was done to real data of double-lined 
spectroscopic binaries, however. The reason is that the individual 
composite observed spectra of all three systems (KIC 11285625, 
KIC 6352430, and KIC 4931738) considered in our work have 
too low S/N for the method to be robust. 

The gssp-Composite software module is much more time con¬ 
suming than the two previously described modules. The reason 
is quite simple: assuming a set up that includes calculation of 
~500 synthetic spectra, a 1 000 A wide wavelength range, 8 
CPUs, and three grid points in RV for each of the components, 
the total number of the grid combinations to be processed by the 
code is close to 1 500 000. This is a factor of almost 10 more than 
for the same set up in the gssp_binary module, which transforms 
into about an hour of computation time. 


3. Summary 

In this paper, we have presented a Fortran 90 code for spectrum 
analysis of single stars and double-lined spectroscopic binary 
systems. The code makes use of the OpenMPI distribution and 
runs in parallel mode. 

There are three independent software modules included, 
each one is designed for the application to a certain type of 
data. The gssp_single module can be applied to the spectroscopic 
data of single stars as well as to the disentangled, light diluted 
spectra of spectroscopic binary components. In the latter case, 
the spectra of the stars are interpreted independently of each 
other and the code assumes wavelength-independent light dilu¬ 
tion for each of the stars. The two other modules, gssp_binary 
and gssp_composite, are specifically designed to fit the spectra 
of double-lined spectroscopic binary systems. Where the first 
module analyzes disentangled spectra, the second one focuses 
on direct fitting of the composite spectra. The gssp_binary mod¬ 
ule fits disentangled spectra of both binary components simulta¬ 
neously and assumes that the light dilution factor of each of the 
components is a function of wavelength. The light factors are re¬ 
placed in the grid by the ratio of the components’ radii, which is 
obviously the same for both stars. This way, an additional con¬ 
straint is put on individual light factors, namely their sum should 
equal unity at each considered wavelength. The algorithm im¬ 
plemented in the gssp_composite module is essentially the same, 
except that composite spectra are fitted instead of disentangled 
ones. Radial velocities of the individual stellar components can 
also be included into the fit if necessary. 

No matter which module is used for the spectroscopic anal¬ 
ysis, the user is strongly advised to start from reasonably coarse 
grid(s) which would cover large parameter ranges. This step is 
absolutely necessary to finally find a global minimum solution 
but not to be stuck in a local minimum. For example, we advise 
to stick to the range in the effective temperature that would cover 
the entire range of spectral type in question. The full ranges for 
the surface gravity and the metallicity currently available in the 
LLmodels library of atmosphere models (see Table|T|for details) 
would be a good choice for the initial run. However, the step 
sizes should be set to rather large values (e.g., 0.5 and 0.4 dex 
for logg and [M/H, respectively]), to keep the calculation time 
at a reasonable level. As soon as the first, coarse grid solution is 
obtained, the user is advises to gradually narrow down the pa¬ 
rameter ranges, also making the grid finer. To speed up the cal¬ 
culations with the gssp_composite module, we would recommend 
to make a priory estimate of the radial velocities of the individ¬ 
ual binary components (if at all possible) and to keep them fixed 
throughout the analysis. 

The gssp_single and gssp_binary software modules are quite 
fast and the calculation time is reasonable even when the code is 
run on a quad core desktop or laptop. If the disentangled spec¬ 
tra of both binary components are available, and the difference 
in spectral types of the two binary components is not large, we 
advise to start the analysis with the unconstrained fitting mode 
implemented in the gssp_single module. This way, very good ini¬ 
tial guesses can be obtained for both stellar components of a 
binary system within short time. With a good initial grid setup 
for the gssp_binary module, the physically more correct solu¬ 
tion for both binary components can be obtained within a lim¬ 
ited amount of time. The disentangled spectra of the compo¬ 
nents of a binary systems can be obtained in different ways; the 
m ost widely used method ologies and codes are those presented 
by ISimon & Sturm (11994 sp ectral disen t anglin g in wavelength 
domain i. lHadraval ( 1995h and llliiic et al.l d200 4. Fourier spectral 
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disentangling). The gssp_composite module is more time con¬ 
suming than the two previous ones, thus we advise to use it when 
i) very good quality data (S/N>150) are available, ii) good ini¬ 
tial guesses for the fitted parameters are provided, and iii) at least 
some of the parameters can be fixed, e.g., from the light curve 
solution. 

Test runs of all three modules on the simulated data have 
shown that all input parameters are well recovered in the ma¬ 
jority of cases. For the fitting of disentangled spectra, S/N>100 
is likely to be sufficient for decent analysis, whereas the com¬ 
posite spectrum fitting generally requires a S/N above 150. The 
gssp software package is also suitable for the analysis of pecu¬ 
liar stars, though one has to be aware of the fact that the global 
metallicity of the stars can be significantly under/overestimated 
should the peculiarity concern one of the dominant chemical el¬ 
ements in the spectrum. 

The application to real data has shown that for binary sys¬ 
tems consisting of two stars with significantly different spectral 
types, the unconstrained fitting is likely to deliver not entirely 
correct parameters. This is due to significant variations of the 
light dilution factor with wavelength, as was the case for the 
KIC 6352430 system analyzed in this study. 

In future, we plan to extend the software to the analysis of hot 
stellar objects (early O-B spectral types) by including non-LTE 
effects into the calculation of synthetic spectra. We will use the 
so-called hybrid approach assuming non-LTE radiative transfer 
calculations for chemical elements of interest, based on classi¬ 
cal LTE atmosphere models. The justification of such approach 
is discussed in lNieva & PrzvbiflaH2007() . The implementation is 
rather straightforward and consists of using pre-computed de¬ 
parture coefficients from LTE in the SynthV radiative transfer 
code. The coefficients can be computed based on the LTE at¬ 
mosphere models with software such as the tlust\0 (iHubenvl 
ll988l l code. We plan to perform such calculations for the en¬ 
tire hot part of our grid of LLmodels atmosphere models, so that 
non-LTE spectrum synthesis can be done for any relevant model 
when it is needed. Computing departure coefficients (including 
elements with complex atomic structure) for a rather large sub¬ 
grid of atmosphere models is a quite time-consuming procedure, 
which may take up to several months of computation time. This 
prevents us from including the non-LTE option in the current 
version of the software release. 
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Table A.l. Results of the application of the gssp_single software module to the simulated data of KIC 11285625 (model 1) and 
KIC 6352430 (model 2) assuming different S/N ratio and dilution factors. The ’’Parameters of the input spectrum 1 ' refer to the exact 
parameter values assumed in the models. The error bars represent lcr uncertainties computed from y 2 statistics and taking into 
account any possible correlations between different parameters. The parameters are given in the form primary/secondary. 


model 

T e ff 

(K) 

logg 

(dex) 

[M/H] 

(dex) 

(kms 

v sin i 
(kms -1 ) 

Dilution 

factor 

S/N 




Parameters of the input spectrum 



60 

150 

la 

7 000/7 200 

4.0/4.2 

-0.2/+0.0 

1.0/2.5 

15/9 

0.95/0.05 

lb 

7 000/7 200 

4.0/4.2 

-0.2/+0.0 

1.0/2.5 

15/9 

0.60/0.40 

60 

150 

lc 

7 000/7 200 

4.0/4.2 

-0.2/+0.0 

1.0/2.5 

15/9 

0.20/0.80 

60 

150 

2a 

12 750/6 800 

4.0/4.3 

-0.1/-0.3 

2.0/3.0 

70/10 

0.95/0.05 

60 

150 

2b 

12 750/6 800 

4.0/4.3 

-0.1/-0.3 

2.0/3.0 

70/10 

0.60/0.40 

60 

150 

2c 

12 750/6 800 

4.0/4.3 

-0.1/-0.3 

2.0/3.0 

70/10 

0.20/0.80 

60 

150 




Results 





la 

7 020/™/7 315/}°° 

A nS+0-08 !A 9C+0.55 

4 -Ho_0.1o/ 4 *^-0.60 

-0.15Z°$/+0.13/°:£ 

1.12/° : }i/2.77/° : 8° 

15.2/«/8.7/|5 

Q 93+0.03/q 05 + 0.01 
- 0 . 03 / U-VA - 0.01 

60 


7 012/}°/7 200Z/°° 

9 00+0.05/4 14 + 0.35 
-0.05'-0.35 

-0.19/° : >0.00/° : |5 

1.02/° : |°/2.53/° : | 

15.1/°;2/9.0/|| 

0 QC+0.02/0 05 \ +0.005 
Kj ‘ yj -0.02/ wuJ1 -0.005 

150 

lb 

7 020/“/7 215/“ 

4 00 +014 /4 13 +0 - 15 

-0.14/^- -0.15 

-0.16/°;°}/ + 0.06/°;°} 

1 03 +017 /2 65 +018 

1 -0.17' z '* u -0.18 

15.2/°;}/9.1/°;} 

0 60 +0 02 /0 41 +0 02 
u. uu _0 02 / 1 - 0.02 

60 


6 995/40/7 195/^f 

9 00+0.06 /4 1 q+0.07 
J,yy -0W^' l - y -0 .07 

-0.22/° : °}/+0.00/° : ^ 

1 A9+0.05 /9 c 9 +0.07 
i,uz, -0.05 /Z - J -0.07 

15 1 +0 - 2 /9 l +0 - 2 
a ^‘ a -0.2/' 7 ‘ a -0.2 

0 61 +001 /0 41 +0 01 
u * UA -o.oi /u, ^ A -o.oi 

150 

lc 

7 070/|^/7 215/5° 

3.85/°;}°/4.16/°;°} 

-0.12/°;!5/ + 0.01/°;°} 

1 IQ+0.40/9 C 7 +O.IO 

A - 1 -0.40/-0.10 

15.2/JI/8.9/°;2 

0.20/°;“/0.79/°;°} 

60 


7 OIO/Oq/7 200/j° 

3.98/°;}}/4.20/°;« 

-0.16/° : °}/ + 0.00/° : °} 

1 04 +016 /2 51 +0 04 

1 -0.16' J -0.04 

15 l +a5 /9 0 +0 - 2 

AJ - - 0 . 5 '^ - 0.2 

0.21/°:°;/0.80/° : °! 

150 

2a 

12710/1^/6 870*™ 

4 O3 +0 09 /4 37 +0 - 58 

- 0 . 08 /°;!J/- 0 . 15 /°;}o 

2.36/°;®/3.36/};™ 

70.0/°;°/9.8/|| 

0 96 +0 03 /o 06 +0 02 

u • ^ u -0.03/ u - uu -0.02 

60 


12 750/“/6790/|™ 

4 on +0 - 04 /4 ?9 + 0 - 46 
^• uu -0.04/^-^-0.50 

-0 10 +0 05 /-0 30 +004 

U,A -0.05' '■'* J -0.04 

2.03/° : g/3.03/°;°} 

70.1/|;|/10. 1/qI 

0 95 +0 01 /0 05 +0 01 
-o.oi /u,u - 0.01 

150 

2 b 

12 860/™/6 810/“ 

4 05 +013 /4 20 +017 

-O-lO^*^ -0.17 

-0 05 +016 /-0 25 +0,08 

u.uj_oi 9 / u.z,j_ 008 

2.55/° : g/3.15/°;l 

69.5/°j/10.1/°;} 

0 60 +0 03 /0 41 +0 02 
u. uu_o 03 / 1 - 0.02 

60 


12 765/|g/6 800/;}° 

4 qq+0.06/4 9Q+0.06 

-0 09 +0 08 /-0 28 +0 - 07 

u * u *-0.08/ u - zo -0.07 

2.00/° : }°/3.04/°;}2 

70. 1 / 40 /IO. I /07 

0 61 +001 /0 40 +0 01 
U - UA -0.01/ U -^ W -0.01 

150 

2 c 

12 860/*°°/6 810/“ 

4.05/°;}°/4.26/°;°} 

-0 01 +a50 /-0 30 +006 

u,vja -0.50 / -0.06 

1.80/} : 5°/3.08/° : ;j 

c+15.0/0 q+0.4 
/w -15.0'^ ,: -0.4 

0.19/°;“/0.81/°;“ 

60 


12 800/§°/6 800/}° 

4.01/° : }5/4.28/° : °} 

-0 10 + °- 22 /-0 30 +0 04 

U - A -0.27' -0.04 

1 qi +1.70/O A9+0.06 

A A —1.70/ -0.06 

70 2 +10 0 /10 0 +0 - 2 
/u *^-10.0 ' aw - u -0.2 

0.20/° : °|/0.80/°;°! 

150 
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Table A.2. Same as Table lA.fl but for the results of the application of the gssp_binary software module. 


model 

Teff 

(K) 

l°gg 

(dex) 

[M/H] 

(dex) 

£ 

(kms ') 

v sin i 
(km s -1 ) 

R 1 /R 2 

S/N 



Parameters of the input spectrum 



60 

150 

la 

7 000/7 200 

4.0/4.2 

-0.2/+0.0 

1.0/2.5 

15/9 

1.45 

lb 

7 000/7 200 

4.0/4.2 

-0.2/+0.0 

1.0/2.5 

15/9 

2.50 

60 

150 

lc 

7 000/7 200 

4.0/4.2 

-0.2/+0.0 

1.0/2.5 

15/9 

4.00 

60 

150 

2a 

12 750/6 800 

4.0/4.3 

—0.1/—0.3 

2.0/3.0 

70/10 

0.50 

60 

150 

2b 

12 750/6 800 

4.0/4.3 

-0.1/-0.3 

2.0/3.0 

70/10 

0.90 

60 

150 

2c 

12 750/6 800 

4.0/4.3 

-0.1/-0.3 

2.0/3.0 

70/10 

1.40 

60 

150 

la 

7 OlO/jjjj/7 280/j3 0 

3 95 +0 - 19 /4 30 +0 - 22 
-0.20' -J -0.26 

Results 

-0.16/0;OS/+0.09/o : “ 

1 05 +0 - 21 /2 45 +0 - 22 

-0.20/^-0.23 

15 2 +0 5 /9 l +0 - 6 

AJ -0.5/^-0.6 

I 47+0.11 

II -0.11 

60 


7 000/“/7 200/“ 

4 00 +0 06 /4 20 +0 08 

-0-20/o : >0.01/r 

0 98+0.06/2 5o+o- 09 
yj ' yo -Q.06 l ^" } -0.08 

15 l +03 /9 0 +0 - 3 

1J -0.3' 7- -0.3 

1 45+0.02 
A- ^ -0.02 

150 

lb 

7 020/“/7 320/^8 

4.00/0^/4.28/0;“ 

-0 16 +007 /+0 16 +018 

V. 1 U_ Q Q7 / TU. 1 U_q 2Q 

I A9+0.17/9 cq+0.46 
A,u -0.17' ^ -0.43 

15.2/04/8.8/11 

9 79+O.25 
^' -0.25 

60 


7 000/^/7 200/“ 

4.00/o : oo/4.20/o;5 

-0 20 +0 03 /+0 00 +006 
-0.03/ -0.06 

i a 1+0.06/9 59+O.I8 
1,Ui -0.06 /Z " J -0.18 

15 2 +0 - 3 /9 0 +0 - 5 

AJ -0.3' 7 -0.5 

9 59+0.14 
-0.14 

150 

lc 

7 020/“/7 320/™ 

4.00/0^/4.28/0“ 

-0 16 +007 /+0 16 +018 

KJ ‘ 1 -0.07' 1 -0.20 

1 m +0 - 17 /? SQ +0 - 46 
1 ,wj -0.17/^ ,J7 -0.43 

15.2/04/8.8/JI 

9 79+O.25 
^ -0.25 

60 


7 015 +60 /7 415 +46 ° 

' u -60/ ' ^ -450 

3.98/o : !!/4.45/| : ™ 

-0.16/0;“/ + 0.15/030 

1 05 +014 /3 30 +1 - 55 

A-U -0.16' DmD -1.46 

15.2/^/8. l/« 

4 Q4+0.8O 
^•^-1.14 

150 

2a 

12590/^0/6810/™ 

4 01 +0.09/4 27+ 017 
^■ ul -0.09 / ^- z,/ -0.17 

-0.08/0;22/-0.28/0;00 

9 99+1.46/9 nc+0.25 
-1.43'-0.25 

72-2/J| ; o/10.2/Og 

0 51+ 0 05 
u,J1 -0.05 

60 


12 750/1^/6 800/™ 

4 OO +0 - 05 /4 9Q+ 0 - 09 
-0.05'^-^ -0.09 

-0 12 +0 09 /— 0 30 +0 ’ 04 

'-'• lz, -o.o9 / -0.04 

1 9 c+0.60/9 00 +01 ° 

1 -^-0.60/-0.10 

71.0/40/10.1/0;! 

0 51 +0 02 
^•^-0.02 

150 

2b 

12 710/|5/6 880/252 

4.04/0;!!/4.58/o : “ 

-0.08/0;|2/-0.27/o : oo 

9 97+O.95/9 99+O.49 

a.j / _ l 27 /j.jj_ 0 60 

70.3/8;9/9.9/[3 

0 96 +0 - 16 

u, ^-0.13 

60 


12 770/^5/6 825/™ 

4.00/0;05/4.28/0 : J2 

-0.12/0;0B/-0.30/0;05 

1 95+0.39/9 o2 +0 - 21 
-0.39' w -0.21 

79 9+3.9/n q+0.6 
-3.9' 7 -0.6 

0 91+0.05 

VJ, ' ? -0.05 

150 

2c 

12 665/205/6 810/^ 

4 O? +0 - 10 /4 1 9+0-85 
^• UZ, -0.10 /H - A -0.85 

-0.08/o : |5/-0.30/o : ™ 

9 9q+I.OO/9 99+I.25 
-1.02' J * JJ -1.00 

70.0/50/9.9/23 

1 40 +0 - 22 

1 -^ w -0.22 

60 


12 775/88/6 790/1™ 

9 qq+ 0.05 /4 94+0.40 
' ? ' :?o -0.05/^' z ' -0.33 

-0.11/o : 02/-0.31/o : |5 

2 00 + °- 39 /3 02 +0 - 45 

z - uu -0.41'- : ’- uz '-0.42 

70.2/3 j/10. i/i ; 2 

1 40+0.O8 

1 -0.08 

150 
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Table A.3. Same as Table lA.fl but for the results of the application of the gssp_composite software module. 


model 

T e tf 

l°gg 

[M/H] 


v sin i 

RV 

R1/R2 

S/N 

(K) 

(dex) 

(dex) 

(kms ') 

(kms -1 ) 

(kms -1 ) 






Parameters of the input spectrum 



60 

150 

la 

7 000/7 200 

4.0/4.2 

- 0 . 2 /+ 0.0 

1.0/2.5 

15/9 

-40/50 

1.45 

lb 

7 000/7 200 

4.0/4.2 

- 0 . 2 /+ 0.0 

1.0/2.5 

15/9 

-40/50 

2.50 

60 

150 

lc 

7 000/7 200 

4.0/4.2 

- 0 . 2 /+ 0.0 

1.0/2.5 

15/9 

-40/50 

4.00 

60 

150 

2 a 

12 750/6 800 

4.0/4.3 

-0.1/-0.3 

2.0/3.0 

70/10 

-40/50 

0.50 

60 

150 

2 b 

12 750/6 800 

4.0/4.3 

—0.1/—0.3 

2.0/3.0 

70/10 

-40/50 

0.90 

60 

150 

2 c 

12 750/6 800 

4.0/4.3 

-0.1/-0.3 

2.0/3.0 

70/10 

-40/50 

1.40 

60 

150 




Results 





la 

7 025 +100 /7 175 +100 

4 02 +a21 /4 11 +0 - 22 
- 0 . 21 '^ ,A a - 0.22 

-0 12 +010 /+0 10 +012 
u,A ^-o.io'™- AVi -o.i 2 

1 12 +a25 /2 60 +0 - 27 
• 1 -0.26'^ ,U -0.27 

1 c 9+O.6 /q 9+O.7 
A ~ , ‘ -0.6' ^ —0.7 

^10.00/°;|/+50.00/°;;; 

1 98+0.10 
^ - 0.10 

60 


6980/®/7210/® 

9 qq+0.06 ia 1 7+O.O8 
' J ‘ yy -0.W^- 1 7 -0.08 

-0 19 +0 03 /+0 02 +004 

1 ^-0.03/ ™ >uz, -0.04 

1 05 +0 09 /2 50 +0 09 

-0.09'-0.09 

15 l +0 - 3 /9 0 +0 - 2 

AJ - l -0.3 /y - yj -0.2 

^40.00/ojl/+50.00/° ; ;4 

1 4 c+0.02 
1 -0.02 

150 

lb 

— 

— 

— 

— 

— 

— 

— 

60 


6 990/4g/7 180/® 

9 Q7+0.06/4 i/:+0.16 

1 -0.16 

-0 20 +0 02 /+0 01 +0 05 
u - z,u -0.02/ ™* u 1 -0.05 

1 00 +0 06 /2 66 +018 
A - UU -0.06 /Z,,U -0.18 

15.111/8-9/gJ 

-39.95/° : |°/+50.07/j; : ;« 

9 r 9+0.06 
-0.06 

150 

lc 

— 

— 

— 

— 

— 

— 

— 

60 


6975/®/7 190/5® 

0 Q7+0.05/4 19+O.4O 
J ‘ y ' -0.40 

-0 21 +0 03 /+0 04 +013 

1 -0.03/ -0.13 

1 01 +0 05 /2 75 + 0 - 50 

i,WA -0.05 / ^- /J -0.50 

i5i/8;l/8.4/;| 

-40.00/q J J/+50.00/q|q 

4 1 1 + 0.20 
^ ,A - 0.20 

150 

2 a 

12990/;®/6850/!® 

4.02/“ : ^/4.27®;| 

-0 22 +0 - 25 /-0 30 +0 07 

-0.25 ' U “ ;,U -0.07 

9 77 +1 - 75 /3 04 +0 - 32 
^ 7 -1.71'-0.32 

73 6 +1L5 /9 8 +0 - 6 

7 ^ -11.5' 7 -0.6 

-40.0/«/+50.00/“ : |; 

0 51 +0 04 

u>ja -0.04 

60 


12 760/JJ°/6 800/|| 

4.00/^/4.28/g;°; 

-0 12 +0 08 /-0 27 +0 05 

W,A -0.08 ' 7 -0.05 

9 nQ+ 0 -58 n m +0 °8 
-o.ss'-^-o.os 

70.6/^|/9.8/o| 

-40.05/; : ®/+49.98/“ : ;; 

0 51+ 0 03 
w,ja -0.03 

150 

2 b 

— 

— 

— 

— 

— 

— 

— 

60 


12 775/| 1 °/6 810/” 

9 Q7+0.05/4 90+0.16 

-0 16 +0 08 /-0 30 +005 

w,au -0.08' -0.05 

1 75+0.50/3 Q4+0.20 
A * 7 -0.50'-0.20 

72.2/«/9.6/“ : ; 

-40.03/|®/+50.00/“;^ 

0 91+0.03 
w> ^ a -0.03 

150 

2 c 

— 

— 

— 

— 

— 

— 

— 

60 


12 775/^/6 790/j® 

9 qo+0.05 M 94+0.40 
• J -^°-o.05 / ^- z ' -0.33 

-0 11 +° 07 /— 0 31 +015 

1 -0.07/ u - :,i -0.15 

2 00 +0 - 39 /3 02 +0 - 45 

z,.uu_04i/^.uz,_042 

70.2/;j/10.1/|“ 

-40.00/S;J|/+50.00/®;!J 

1 40+0.08 

1 -0.08 

150 
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km s ' 1 

0 96 +aos 
yl - yu -om 

O.S2/« : ™ 

2 01 +al9 
z '- ul - 0.21 

9 4C+0.52 
^*^-0.55 

0 96+ 0 07 

u.?u 007 

0.63l;l 

9 qo+ 0.50 
— u -0.50 

2 00 +0 - 50 
-0.50 

1 01 +0 ° 5 

1 ■ ul -0.05 

0 97 + 0.08 
^ y1 -0.08 

2 07+0.42 
^ ' -0.42 

2 Q 7 + 0.42 
^ ' -0.42 

v sin i 

km s _1 

15 l +a2 

15-1^1 

70 3 +4 -° 

'0-9_4.q 

69 5 +4 0 

-4.0 

15 l +a2 

IJ.I -0 2 

15 . 0 I 3 

70 1 +43 

' 0 - 1 -4.0 

70 1 +43 

/u ‘ a -4.0 

15 l +a2 
IJ. 1-0 2 

15 l +0 - 2 

1J ‘ 1 - 0.2 

70.6/ 3 :* 

70.6/ 3 :® 

[M/H] 

dex 

- 0 . 2 * 

-0 05 +0 07 

0-0 J -0.07 

- 0 . 1 * 

-0 15 +009 

u - 1 -0.09 

- 0 . 2 * 

0 02 +012 
- 0.12 

- 0 . 1 * 

-0 12 +012 
iZ, - 0.12 

- 0 . 2 * 

A AC+0.06 

U.UD_o 06 

- 0 . 1 * 

-0 ll +0 07 
u - 1 -0.07 

[El/H] 

dex 

_4 48+ 0 < >5 
^• ho -0.05 

a 45 + 0.08 

-1 30 + ° 06 

1 " ; ’O-0.07 

1 ^ 1 +0-07 

1 • ^ 1 -0.07 

_4 48 + 0.03 
^• ho -0.03 

a 49 + 0.09 
-0.09 

1 9 Q+O.O 7 

A -^-o.io 

1 9Q+0.07 
L -^ y -o.io 

_4 48+ 0 02 
^• ho - 0.02 

A 4 ^+0.05 
-0.05 

1 9 A+O.O 6 
L ’^ y -0.01 

1 9Q+0.06 

1 ' 46 -0.07 

R 1 /R 2 


9 c 9+0.05 
-0.05 

9 45+ 042 
— H -0.12 

0 gQ+ 0 . 0 4 
u - y O_ 004 

0 90 + ° 05 
''■^-0.05 

n ^n+o.oi 

0.60_ooi 

n c/^+0.07 
-0.07 

0 61 +0 01 
u ’ ui - 0.01 

0 61 +0 01 
u ’ ul - 0.01 

9 54+0.06 
-0.06 

9 oc+0.08 
z, - z,J -0.08 

0 9l +0 04 
u ‘^ a -0.04 

0 9l +0 04 
u> ^ a -0.04 

RV 

km s _1 

— 

— 

— 

— 

— 

— 

— 

— 

-40 O 94009 

— _ 0 09 

-39 98 +010 

jy. y o_ 0lQ 

-40.06/ 2 :® 

-40.03/1 







Secondary 







T c s 

K 

7190/™ 

7 185/1“ 

6 810/1 

6 800/1 

7 195 +40 

7 195 +40 

' 1:7 -40 

— 

— 

7185/g 

7190/1 

6 805/1 

6 805/1 

log g 

dex 

4 91+0.14 
^*^ A -0.16 

+ 1611 

4 28 +017 

^• z,o -0.17 

4 99+0.16 
-0.18 

4 1Q+0.06 

1 -0.06 

4 iQ+0.06 

1 -0.06 

— 

— 

4 1 5+0.16 
-0.16 

4 35 +017 

-0.17 

4 9C+O.I6 
^•^-0.18 

4 95+0.16 
-0.18 


km s _1 

9 C9+O.I7 

0 ^Q+0.24 
^ -0.39 

3 02 +0 - 22 
-0.20 

2 96 +a21 

z " ,u -0.25 

9 c 9+0.06 
-0.06 

9 c 9+0.06 
-0.06 

— 

— 

2 64+0-22 
z " u -0.22 

2 05 +a25 

0.24 

3.071™ 

3.0?1™ 

v sin i 

km s _1 

9 o +0 - 5 

y -' J - 0.5 

0 4+0.7 
'*■ -0.7 

101 +0 - 5 

AU,A -0.5 

10 l +a5 

lu --0.5 

9 1+0.2 
7-0.2 

9 1+0.2 
^-0.2 

— 

— 

8 9 +03 

°- y - 0.5 

A 4+0.7 
-0.7 

10 . 01 J 

10 . 0 IJ 

[M/H] 

dex 

0.0* 

-0 05 +0 09 
O' 1 -' -0.09 

a 9 0+0.08 
-0.08 

-0 30 +al0 
-0.10 

0.0* 

0 oo +003 

-0.03 

— 

— 

0.0* 

-0-0511 

a 90 +0.08 
U - Z °-0.08 

a 90+0.08 
U - Z8 -0.08 

[El/H] 

dex 

-5 01 +035 

-5 05 +a52 

-0.57 

— 

— 

-4 98 +014 
^• 7 o - 0.18 

a 98+0- 14 
^■ ?o -0.18 

— 

— 

-5 05 4034 

J.UJ-034 

-5-0611 

— 

— 

R1/R2 


9 59+0.05 
^ -0.05 

9 45+0- 12 
• J -0.12 

0 9o+° 04 
u -* u -0.04 

0 90 + ° 05 
u,?u -0.05 

0 40+ 0 01 
001 

0 40+ 0 01 
u.hu 001 

— 

— 

9 54+O.O6 
• J -0.06 

9 oc+0.08 
^ -0.08 

0 91 4004 

u - yl -0.04 

0-911S 

RV 

km s _1 

— 

— 

— 

— 

— 

— 

— 

— 

49 98+ 011 
^ y - yo -0.11 

50-0011 

50.001 24 

50.05 lU 
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A. Tkachenko-: Grid Search in Stellar Parameters: a software for spectrum analysis of single stars and binary systems 

Appendix B: The gssp package: Installation, running, and configurations files 

B. 1. Installation and running 

The Grid Search in Stellar Parameters (gssp) package is a compilation of three separate modules, each one designed for the analysis 
of a certain type of spectroscopic data. The code makes use of an open source Message Passing Interface (OpenMPI) implementation, 
which has to be installed along with the Intel Fortran Compiler on the machine where the gssp code is supposed to run. The version 
of OpenMPI 1.3.3 from July 14th, 2009 is suitable; any stable version of the Intel Compiler is good but the release should not be 
older than the one from the 12th of February, 2012 (version 12.1.3). 

The code compiles with the following simple command: 

mpif9® -o executable_name *.f* 

where“executable_name” is the name of an executable file (can be anything the user likes). The compilation command is the same 
for all the package modules. For the compilation to be successful, one needs to run the compilation command twice. During the 
first run, the GSSP_module.f90 file is processed and the *.mod file with the corresponding name is produced. This file contains the 
majority of declarations as well as some other relevant information for other subroutines. During the second run of the compilation 
command, the compiler gets the information from the *.mod file and all remaining Fortran source files get successfully processed. 
To run the code, one need to use the following command: 

mpirun -n N [-hostfile Host_file_name] executable_name Input_file_name 

where N is the number of CPUs. The arguments given in brackets are optional and not needed when the code is run on a single 
machine with multiple-cores. The arguments become useful when the code runs on a cluster PC with several nodes involved. The 
host file contains information about the names of the individual nodes on which the code will be running, and the number of CPUs 
on each of those nodes. An example of such hostfile is given below. 

nodeOl slots=8 
node02 slots=8 
node03 slots=8 
node04 slots=8 
node05 slots=8 


B.2. Configuration files 

In this section, we discuss the structure of the input files for all three modules of the gssp package. Our intention is to give as much 
detail as needed for each entry from those files. 

The gssp_single module 

Teffstart Tejfstep Teff-end ! effective temperature range (K) 

loggstart loggstep logg-end ! log surface gravity range (dex) 

vmicrostart vmicrostep vmicrojend ! microturbulent velocity range (km s _1 ) 

vsinistart vsinistep vsinijend ! projected rotational velocity range (km s _1 ) 

dilution-flag factorstart factorstep factorsnd ! light dilution flag and light dilution factor range 

abundflag [M/H]start [M/H]step [M/H]_end ! abundance flag and metallicity range (dex) 

element-id abundstart abundstep abundsnd ! chemical element id (e.g., Fe or He) and abundance range 

vmacro resolution ! macroturbulent velocity (km s 1 ) and resolving power 

abundance .path ! absolute path to the abundance table(s) 

model-path ! absolute path to the library of atmosphere models 

vmicrojnodel mass-model ! atmosphere model microturbulent velocity (kms -1 ) and mass (M 0 ) 
modeLcompflag ! atmosphere model chemical composition flag (ST, CNm, or CNh) 
nranges wave step mode ! Number of wavelength regions, wavelength step (A), operational mode flag 
spectrum-path ! Path to the observed spectrum 

RV-factor contin-factor RV RVflag ! RV scaling factor, continuum cutoff factor, RV value, RV option 
wave start wavesnd ! wavelength range(s) in A 
wavestart wave-end 

The first four entries in the configuration file set up the ranges and step sizes for the effective temperature (7’ c iij, logarithm of 
surface gravity (log g), microturbulent velocity (f), and projected rotational velocity (v sin i). For example, 
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6800 100 7200 ! effective temperature range (K) 

2.9 0.2 3.9 ! log surface gravity range fdex) 

2.8 0.4 4.4 ! microturbulent velocity range (km s -1 ) 

35 5 50 ! projected rotational velocity range (km s~') 

The next two entries set up the grid ranges for the light dilution factor and the global metallicity [M/H] of the star. In both cases, 
certain flags precede the parameter ranges. Should the light dilution option be used in the calculations (in case of the analysis of 
disentangled spectra of binary components), the value “adjust” must be set for the dilution_flag parameter. Otherwise (analysis 
of single stars), the value of the dilution_flag parameter needs to be set to “skip”. Similarly, the abund_flag takes one of the 
two following values: “skip” or “adjust”. The “skip” option assumes that the global metallicity of the star can be optimized along 
with the other fundamental parameters (7j.fr, log g, f, and v sin /). In this case, individual abundances of all metals will be scaled by 
the same amount. The “adjust” option assumes that the metallicity parameter will be replaced in a grid by a chemical element (any, 
except hydrogen) abundance. The individual abundances are optimized for one chemical element at the time; the metallicity value 
[M/H] must be fixed while all other fundamental parameters can be kept free. For example, 

adjust 0.6 0.1 1.0 ! light dilution flag and light dilution factor range 
skip -0.1 0.1 0.2 ! abundance flag and metallicity range (dex) 

For detailed chemical composition analysis of the star, we recommend to do the analysis in two steps: 1) optimizing the global 
metallicity value together with the other four (and, optionally, the light dilution factor) fundamental parameters; 2) fixing the 
metallicity to the value derived in the first step and optimizing individual abundances along with the other parameters. 

The entry element Jd abundstart abundstep abuiuLend is only relevant when the individual abundances have to be optimized. 
Below are the entry examples for Fe: 

Fe -4.59 0.05 -4.39 ! chemical element id (e.g., Fe or He) and abundance range 
and for He 

He 0.0783 0.0005 0.0813 ! chemical element id (e.g., Fe or He) and abundance range 

The element_id parameter refers to the element designation (e.g., Fe for iron. He for helium); the numbers give the range and 
a step width. Should the value of the abund_flag parameter be set to “skip”, the entire entry will be ignored by the code. Note: 
though abundances of all metals are in logarithmic scale, those of hydrogen and helium are not. Thus, when optimizing the helium 
abundance, one has to remember that these are positive values and the range should always be given from the smallest to the largest 
value (see example above). 

The vmacro resolution entry does not need any specific remarks as it gives the values of macroturbulent velocity (in km s') and 
of the resolving power of the instrument. Both parameters are used by the convolution code to add the corresponding broadening to 
synthetic spectra. For example, 

0.0 32000 ! macro turbulent velocity (km s -1 ) and resolving power 

The next two entries, abundance .path and model-jtath , specify absolute paths to a file with individual abundances and to the 
folder containing atmosphere models, respectively. For example: 

/home/UserName/Abundance_table.abn ! absolute path to the abundance table(s) 

/home/UserName/LLmodels/ ! absolute path to the library of atmosphere models 

or 

/home/UserName/abundances/ ! absolute path to the abundance table(s) 

/home/UserName/LLmodels/ ! absolute path to the library of atmosphere models 

A library of LLmodels atmosphere models is included into the the gssp package; the library of interpolated Kurucz models with the 
resolution matching the one of the LLmodels grid is provided upon request. The file with the individual abundances can have any 
name but the extension *.abn is fixed. The program will not run in the abund_f lag = ad j ust mode if the file name with individual 
abundances is not provided. At the same time, the global metallicity parameter [M/H] must be fixed when the name of *.abn file is 
provided in the configuration file. Whenever [M/H] is set as a free parameter, the path should be provided to the folder that contains 
abundance tables corresponding to different global metallicity values (e.g., /home/UserName/abundances/). These global metallicity 
abundance tables have been pre-computed and are included in the gssp package distribution. Summarizing, the abundance_path 
parameter is closely connected to the [M/H] parameter (whether it is free or not) and to the abund_flag: 

- [M/H] is a free parameter: the absolute path to the folder with pre-computed chemical composition tables should be provided 
(e.g., /home/UserName/abundances/). Providing the name of a *.abn file in this case will result in an error message. 

- [M/H] is fixed and abund_flag = skip: the other four fundamental parameters (7/.n, logg, </, and v sin i) and (optionally) the 
dilution factor can be optimized based on the provided table with individual abundances. There are two options here: either one 
provides a path to the file that contains specific abundances (e.g., /home/UserName/Abundance_table.abn), or a path to the folder 
containing pre-computed abundance tables corresponding to different metallicity values (e.g., /home/UserName/abundances/). 
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- [M/H] is fixed and abund_flag = adjust: the other four fundamental parameters (T e ff, log#, and v sin i) and (optionally) 
the dilution factor can be optimized simultaneously with the individual abundances (one element at the time). In this case, a 
path to the file that contains specific elemental abundances has to be provided (e.g., /home/UserName/Abundance_table.abn). 
Otherwise, the code will give an error message. 

The next two entries in the configuration file, vmicro jnodel mass jnodel and model -comp -flag, set the values of the atmosphere 
model microturbulent velocity (in km s 1 ) and mass (M e ), and the value of the atmosphere model chemical composition flag. For 
example,: 

2 1 ! atmosphere model microturbulent velocity (km s ~ 1 ) and mass (M Q ) 

ST ! atmosphere model chemical composition flag (ST, CNm, or CNh) 

All atmosphere models in the provided LLmodels grid have been computed for the fixed value of the microturbulent velocity of 
2 kms -1 . However, the user is free to use his/her own grid of models (in Kurucz format!) computed for a different value of the 
microturbulence. The other two parameters, mass_model and model_comp_flag, were specifically included should the user want 
to use marcsQ (Gustafsson et al. 2008) atmosphere models computed assuming a different mass and/or chemical composition. 
Should the user choose to use the provided grid of the LLmodels atmosphere models, both parameters must be fixed to the values 
given in the above example. Note: when using spherical marcs models, one has to be aware that a hybrid approach is adopted 
in the spectrum analysis - spherical models and a plane-parallel radiative transfer formalism in the spectrum synthesis code. The 
parameter model_comp_flag takes one of the three following values: “ST” (stands for standard), “CNm” (stands for moderately 
CN-cycled), or “CNh” (stands for heavily CN-cycled). The “standard” mixture reflects the typical elemental abundance ratios in 
stars as a function of metallicity in the solar neighbourhood. Two types of CN-cycled marcs models reflect different carbon isotopic 
ratios ( 12 C/ I3 C; see marcs website for details). 

The nranges wavestep mode entry sets the number of wavelength ranges to be considered in the analysis, the wavelength step 
width (in A), and the operational mode. For example: 

2 0.0376 fit ! Number of wavelength regions, wavelength step (A), and operational mode flag 

The mode parameter can be assigned one of the following key words: “fit” or “grid”. The first, fitting mode, assumes that the 
observed spectrum is provided (see below) and a grid of synthetic spectra will be fitted to it. The quality of the fit is evaluated based 
on the ^-criterion; the x 1 distributions are provided for all free parameters as an output (see below). The number of wavelength 
regions can be set to any integer number, the ranges themselves are provided further in the configuration file. The value of the step 
width in wavelength is ignored in this case, and is computed directly from the observations (see below). The second value that the 
mode parameter can take is “grid”. In this case, no fitting to the observations is performed and a grid of synthetic spectra is computed 
in the required parameter space. This option assumes that the spectra are calculated in a single wavelength range, thus the number 
of wavelength ranges (nranges) should be set to unity. The step width in wavelength is also an important parameter in this case, as 
it defines the wavelength grid for synthetic spectra. 

The spectrum-path entry provides a path (absolute or relative) to the observed spectrum. This entry in the configuration file will 
be ignored should the user run the code in the “grid” mode. For example: 

home/UserName/Obs_spectrum.dat ! Path to the observed spectrum 

The observed spectrum is an important part of the fitting mode, thus an error message will be given when the spectrum is not 
provided but the operational mode is set to “fit”. The observed spectrum should be provided in a two-column ASCII file, where the 
first and the second columns refer to wavelength (in A, linear scale) and normalized flux, respectively. Note: the wavelength scale of 
the observed spectrum should be equidistant. The step width in wavelength that will be used for the calculation of synthetic spectra 
is computed from the observations. 

The RV^factor continjactor RV RV-flag entry is relevant when the fitting operational mode is chosen. Otherwise, it will be 
ignored by the code. For example: 

0.50 0.99 -35.5 fixed ! RV scaling factor, continuum cutoff factor, RV value, RV option 

The gssp_single module has an option to compute the cross-correlation function (CCF) between the observations and the first 
synthetic spectrum from the grid. The CCF will be computed when the RV_flag is assigned the key word “adjust”. Otherwise 
(RV_flag = fixed), the code assumes that the radial velocity of the observed spectrum is known and the RV value preceding 
the RV_flag in the configuration file (-35.5 km s in the example above) is used to correct the observations. Should the CCF 
calculation be requested by the user, the RV of the observed spectrum is determined as the first order moment of the CCF. The 
velocity range in the CCF from which the radial velocity is computed is defined by the user in a parameterized way. The control 
parameter is RV_factor (takes the value of 0.5 in the example above), which sets the intensity cut limit in the normalized to a unity 
CCF function. We strongly recommend to use a value of 0.5 or slightly larger, which in practice means that the CCF will be cut at its 
half maximum intensity and the part above will be used to compute the first order moment. Finally, the contin_factor parameter 
takes a value between 0 and 1 and reflects the cutoff value in the normalized observed flux. All flux points above this value (0.99 in 
the example above) will be taken into account in the calculation of the global continuum correction factor. The information about 
the global continuum position is taken from the synthetic spectra; the factor is computed by means of the least-squares fit of the 

4 http://marcs.astro.uu.se/ 
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observations to the theoretical spectrum. In practice, the value of this parameter should be fairly large (between 0.95 and 1.0) to 
ensure that mainly the continuum points are taken into account and the result is not biased by the inclusion of weak spectral lines. 
Should the user not want to do any continuum correction, the parameter can be set to some unreliably large value (e.g., 10000), 
which will force the correction factor to be identical to unity. 

The wavestcirt wave-end entry refers to the wavelength range to be considered in the analysis. The number of entries should 
be equal to the number of wavelength regions (nranges) specified earlier in the configuration file (see above). For example: 

4750 5000 ! wavelength range(s) in A 
5200 5700 

Only one wavelength range can be provided when the operational mode is set to “grid”. In the fitting mode, the number of wavelength 
regions is not limited, which makes the fitting of individual lines in the spectrum possible. 

The gssp_binary module 

number .components ! number of stellar components 

Teff.start 1 Teffstep 1 Teff .end 1 Teffstart 1 Teffstep 2 Teff .end 2 

loggstart 1 loggstep 1 loggjend 1 loggstart 2 loggstep 2 logg-end 2 

vmicrostart 1 vmicrostep 1 vmicrojend 1 vmicrostart 2 vmicrostep 2 vmicrojend 2 

vsinistart 1 vsinistep 1 vsinijend 1 vsinistart 2 vsinistep 2 vsinijend 2 

abundjiag 1 [M/H]start x [M/H]step x [M/H].end 1 abund-flag 2 [M/H]start 2 [M/H]step 2 [M/H].end 2 

elementJd 1 abundstart' abundstep 1 abund-end 1 elementdd 2 abundstart 2 abundstep 2 abund-end 2 

radii-ratio start radii-ratiostep radii -ratio .end ! ratio of the components’ radii 

vmacro 1 vmacro 2 resolution 

abundance ^rath 1 

abundance 4 >ath 2 

model ^>ath 

vmicrojnodel massjnodel 
model .comp Jlag 

nranges ! Number of wavelength regions to be considered in the analysis 
spectrumjrath 1 ! Path to the disentangled observed spectrum of the primary 

spectrum jrath 2 ! Path to the disentangled observed spectrum of the secondary 
RV .factor 1 RV 1 RVfllag 1 contin .factor' RVffactor 2 RV 2 RV.flag 2 contin.factor 2 
wavestart wave.end 
wavestart wave.end 

Since the structure of the configuration file itself and the meaning of the input parameters are very similar to those of the gssp_single 
module, we only discuss extra parameters specific to the gssp_binary software module. In the above example of the configuration 
file, superscripts “1” and “2” refer to the primary and secondary components of a binary, respectively. 

The number.components refers to the number of stellar components in the system. Currently, the value of this parameter always 
has to be set to 2. 

radii.ratio start radii.ratio step radii.ratio-end is another entry in the configuration file that needs some attention. The entry 
refers to the binary components’ radii ratio and is obviously the same for both stellar components. Depending on the binary system 
considered, this ratio can have values close to one or above ten. 

The gssp_composite module 

number-components ! number of stellar components 

Teff start 1 Teffstep 1 Teff.end 1 Teff start 2 Teffstep 2 Teff.end 2 

loggstart 1 loggstep 1 loggjend 1 loggstart 2 loggstep 2 logg-end 2 

vmicrostart 1 vmicrostep 1 vmicro.end 1 vmicrostart 2 vmicrostep 2 vmicrojend 2 

vsinistart 1 vsinistep 1 vsini.end 1 vsinistart 2 vsinistep 2 vsinijend 2 

abund-flag 1 [M/H]start x [M/H]step' [M/H].end 1 abund-flag 2 [M/H]start 2 [M/H]step 2 [M/H].end 2 

element-id 1 abundstart 1 abundstep 1 abund-end 1 element-id 2 abundstart 2 abundstep 2 abund-end 2 

radii .ratio start radii .ratio step radii.ratio.end ! ratio of the components’ radii 

RVstart 1 RV step 1 RV-end 1 RVstart 2 RV step 1 RV.end 2 ! radial velocity ranges (km s ') 

vmacro 1 vmacro 2 resolution 

abundance j>ath 1 
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abundance .path. 2 
model-path 

vmicrojnodel massjnodel 
model jcompjlag 

manges ! Number of wavelength regions to be considered in the analysis 
spectrum-path ! Path to the composite observed spectrum of the binary system 
con tin factor ! continuum cutoff factor 
wavestart wave-end 
wavestart wave-end 

The meaning of the majority of the entries is the same as in the gssp_single and gssp_binary modules, thus a detailed discussion is 
omitted here. 

RVstart 1 R Vstep 1 RVsnd 1 RVstart 2 RVstep 2 RV-end 2 is the only entry that requires some attention as it is new 
compared to all previously discussed. The entry obviously refers to the individual component grids in radial velocities (in kms -1 ). 
Should the user want to fix the RV of either of the component stars, the RV_start parameter must be set equal to the RV_end 
parameter. 

B.3. Output files 

All output files are stored in the “outputjfles” folder that is located in the local work directory and is created automatically by the 
code. The first important output file has the name “ModelOverview.txt”. Whenever the user is out of the parameter range provided 
by the library of atmosphere models, the code will exit with an error message, referring the user to this file. It gives an overview of 
all models that do (or do not) exist in the specified parameter space. The file has the following structure: 

[M/H]=-0.1 Teff= 6800.0 logg=2.90 Vmicro=2.0 Mass=1.0 OK 
[M/H]=-0.1 TefF= 6800.0 logg=3.10 Vmicro=2.0 Mass=1.0 OK 
[M/H]=-0.1 Teff= 6800.0 logg=3.30 Vmicro=2.0 Mass=1.0 OK 
[M/H]=-0.1 Teff= 6800.0 logg=3.50 Vmicro=2.0 Mass=1.0 OK 
[M/H]=-0.1 Teff= 6800.0 logg=3.70 Vmicro=2.0 Mass=1.0 OK 

[M/H]=-1.0 Teff= 6800.0 logg=2.90 Vmicro=2.0 Mass=1.0 Not found 

In this particular example, the code informs the user that one of the low-metallicity atmosphere models is missing in the library. 
Obviously, one has to reconsider the parameter space to stay within the ranges of fundamental parameters assumed by the library, 
or to consider using his/her own library of models that covers the parameter range of interest. 

Should the user request calculation of the CCF between the observations and theoretical spectrum (gssp_single and gssp_binary 
modules), the CCF itself will be saved in the “output_files/CCF.dat” file. This is a two-column ASCII file, with the columns 1 and 
2 giving radial velocity and the cross-correlation coefficient, respectively. The radial velocity value determined as the first order 
moment of the CCF is displayed on the screen. No file with the CCF is produced by the gssp_composite module as the RVs of both 
binary components are parts of the grid search parameter space. 

Each of the gssp package modules produces output files that contain the observations and the best fit synthetic spectra. These 
are the files that have to be plotted to visually investigate the quality of the fit. The output observed spectrum is corrected for the 
RV shift (gssp_single and gssp_binary modules) and for possible imperfections in global normalization. This is a two-column ASCII 
file, where the first column gives the wavelength in A and the second one gives the normalized flux. The best fit synthetic spectra 
are also stored in ASCII files, and depending on the software module used, can contain up to six columns. Below we give some 
guidance on the content of the synthetic spectra output files: 

- gssp_single module: in the case of a single star spectrum fitting, the output file with the synthetic spectrum will have an extension 
*.rgs. The only relevant data for the user are in the first and the second columns, which give the wavelength in A and normalized 
flux, respectively. Should the user want to fit the disentangled spectrum of a binary component in the unconstrained mode 
(light dilution mode), the final light diluted synthetic spectrum is stored in a two-column ASCII file with *.rgs extension, which 
contains the wavelength in A (column 1) and the normalized diluted flux (column2). 

- gssp_binary module: in the constrained fitting mode, the output best fit synthetic spectra files (one per stellar component) 
contain three columns: 1. the wavelength in A; 2. the normalized synthetic flux; 3. the normalized light diluted synthetic flux. 
Thus, columns 1 vs. 3 have to plotted for direct comparison with the observed disentangled spectra of the individual stellar 
components. 

- gssp_composite module: in the case of the composite spectrum fitting, the best fit synthetic spectrum is stored in a four-column 
ASCII file. The meaning of the columns is the following: 1. the wavelength in A; 2. the composite normalized synthetic flux; 3. 
the normalized synthetic flux of the primary; 4. the normalized synthetic flux of the secondary. For direct comparison with the 
observations, the columns 1 vs. 2 have to be plotted. 
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Depending on the gssp software module used, the output files containing the spectra have slightly different names. In any case, 
the file names speak for themselves and the observations can easily be distinguished from the best fitting synthetic spectra. 

The;p 2 master file is saved in the “output_files/Chi2_table.dat” ASCII file. Depending on the gssp software module used, the file 
can contain up to 19 columns. Below the headers of the x 1 output files from the different gssp modules are given. 

gssp_single module, single star spectrum: 

[M/H] r eff logg £ [abundance] v sin i ^ 2 _inter contin_factor reduced^ 2 x\ <T 

gssp_single module, disentangled spectrum of a binary components: 

[M/H] 7/n log g £ [abundance] v sin i dilution_factor ^f 2 _inter contin_factor reduced^ 2 x\ 

The relevant columns for the user are all fundamental parameters (and, optionally, dilution factor), the reduced^ 2 , and the lcr level 
i 11 x 1 ■ The parameter vs. reduced^ 2 columns can be used to plot individual^ 2 distributions; the x] a level can be used to estimate lcr 
uncertainties. The parameter in brackets is optional and will appear in the final table only if the individual abundances are optimized. 

The x 2 master file is sorted in a way that the best fit solution appears in the first row of the file. The best fit parameters are also 

displayed on the screen. 

gssp_binary module: 

[M/H],- T eSi log gi & v sin i, [abundance,] (radii_ratio) contin_factor, reduced^ 2 0p 2 um x] (r ) 

gssp_composite module: 

[M/H], T e ffj logg/ £ v sin ij [abundance,] RV,- (radii_ratio contin_factor reduced^ 2 x\ o-) 

In the above two cases, the subscript refers to the stellar component (1 - for the primary, 2 - for the secondary). The parameters 
given in parentheses appear in the file only once and always follow the atmospheric parameters of the secondary. In the above exam¬ 
ple of a gssp_binary module output file, the parameters in the^- 2 master file have the following order: 1 ) the atmospheric parameters 
of the primary+continuum correction factor and reduced^ 2 for the corresponding observed spectrum; 2) the atmospheric parameters 
of the secondary, including the components’ radii ratio+continuum correction factor and reduced^ 2 for the corresponding observed 
spectrum; 3) the sum of the individual reduced^ 2 values (x\ +X%) anc * 1°" l eve l in tLih- 

Note: when running the gssp_single software module in the “grid” operational mode, the code produces an additional output 
folder named “rgs_files”. This folder contains the entire grid of the synthetic spectra. 
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